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ABSTRACT 

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 lu- 
minous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72h~^ Gpc"^ over 
3816 square degrees and 0.16 < z < 0.47, making it the best sample yet for the study of large-scale 
structure. We find a well-detected peak in the correlation function at 100h~^ Mpc separation that is an 
excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic 
oscillations on the low-rcdshift clustering of matter. This detection demonstrates the linear growth of 
structure by gravitational instability between z « 1000 and the present and confirms a firm predic- 
tion of the standard cosmological theory. The acoustic peak provides a standard ruler by which we 
can measure the ratio of the distances to z = 0.35 and z = 1089 to 4% fractional accuracy and the 
absolute distance to z = 0.35 to 5% accuracy. From the overall shape of the correlation function, we 
measure the matter density ilmh^ to 8% and find agreement with the value from cosmic microwave back- 
ground (CMB) anisotropics. Independent of the constraints provided by the CMB acoustic scale, we 
find n„i = 0.273 ± 0.025 + 0.123(1 + wq)+ 0.137rJif . Including the CMB acoustic scale, we find that the 
spatial curvature is ^Ik = — 0.010 ± 0.009 if the dark energy is a cosmological constant. More generally, 
our results provide a measurement of cosmological distance, and hence an argument for dark energy, 
based on a geometric method with the same simple physics as the microwave background anisotropics. 
The standard cosmological model convincingly passes these new and robust tests of its fundamental 
properties. 

Subject headings: cosmology: observations — large-scale structure of the universe — distance scale — 
cosmological parameters — cosmic microwave background — galaxies: elliptical and 
lenticular, cD 
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1. INTRODUCTION 

In the last five years, the acoustic peaks in the cosmic 
microwave background (CMB) anisotropy power spectrum 
have emerged as one of the strongest cosmological probes 
(Miher et al. 1999; de Bernardis et al. 2000; Hanany et al. 
2000; Halvcrson ct al. 2001; Lcc ct al. 2001; Ncttcrficld ct 
al. 2002; Bcnoit et al. 2003; Pearson ct al. 2003; Bennett et 
al. 2003). They measure the contents and curvature of the 
universe (Jungman et al. 1996; Knox & Page 2000; Lange 
et al. 2001; Jaffe et al. 2001; Knox et al. 2001; Efstathiou et 
al. 2002; Percival et al. 2002; Spcrgcl ct al. 2003; Tcgmark 
et al. 2004b), demonstrate that the cosmic perturbations 
are generated early {z » 1000) and are dominantly adia- 
batic (Hu & White 1996a,b; Hu et al. 1997; Peiris et al. 
2003; Moodlcy ct al. 2004), and by their mere existence 
largely validate the simple theory used to support their 
interpretation (for reviews, see Hu et al. 1997; Hu & Do- 
delson 2002). 

The acoustic peaks occur because the cosmological per- 
turbations excite sound waves in the relativistic plasma of 
the early universe (Peebles & Yu 1970; Sunyaev & Zel'dovich 
1970; Bond & Efstathiou 1984, 1987; Holtzmann 1989). 
The recombination to a neutral gas at redshift z w 1000 
abruptly decreases the sound speed and effectively ends 
the wave propagation. In the time between the forma- 
tion of the perturbations and the epoch of recombination, 
modes of different wavelength can complete different num- 
bers of oscillation periods. This translates the character- 
istic time into a characteristic length scale and produces a 
harmonic series of maxima and minima in the anisotropy 
power spectrum. 

Because the universe has a significant fraction of baryons, 
cosmological theory predicts that the acoustic oscillations 
in the plasma will also l)e imprinted onto the late-time 
power spectrum of the non-relativistic matter (Peebles & 
Yu 1970; Bond & Efstathiou 1984; Holtzmann 1989; Hu 
& Sugiyama 1996; Eisenstein & Hu 1998). A simple way 
to understand this is to consider that from an initial point 
perturbation common to the dark matter and the baryons, 
the dark matter perturbation grows in place while the 
baryonic perturbation is carried outward in an expanding 
spherical wave (Bashinsky & Bertschinger 2001, 2002). At 
recombination, this shell is roughly 150 Mpc in radius. 
Afterwards, the combined dark matter and baryon per- 
turbation seeds the formation of large-scale structure. Be- 
cause the central perturbation in the dark matter is dom- 
inant compared to the baryonic shell, the acoustic feature 
is manifested as a small single spike in the correlation func- 
tion at 150 Mpc separation. 

The acoustic signatures in the large-scale clustering of 
galaxies yield three more opportunities to test the cos- 
mological paradigm with the carly-univcrse acoustic phe- 
nomenon: 1) it would provide smoking-gun evidence for 
our theory of gravitational clustering, notably the idea 
that large-scale fluctuations grow by linear perturbation 
theory from z ~ 1000 to the present; 2) it would give 
another confirmation of the existence of dark matter at 
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z ^ 1000, since a fully baryonic model produces an effect 
much larger than observed; and 3) it would provide a char- 
acteristic and reasonably sharp length scale that can be 
measured at a wide range of redshifts, thereby determining 
purely by geometry the angular-diametcr-distance-redshift 
relation and the evolution of the Hubble parameter (Eisen- 
stein, Hu, & Tegmark 1998; Eisenstein 2003). The last ap- 
plication can provide precise and robust constraints (Blake 
& Glazebrook 2003; Hu & Haiman 2003; Linder 2003; Seo 
& Eisenstein 2003; Amendola et al. 2004; Dolney, Jain, & 
Takada 2004; Matsubara 2004) on the acceleration of the 
expansion rate of the universe (Riess et al. 1998; Perlmut- 
ter et al. 1999). The nature of the "dark energy" causing 
this acceleration is a complete mystery at present (for a 
review, see Padmanabhan 2004), but sorting between the 
various exotic explanations will require superbly accurate 
data. The acoustic peak method could provide a geomet- 
ric complement to the usual luminosity-distance methods 
such as those based on type la supernovae (e.g. Riess et 
al. 1998; Perlmutter et al. 1999; Knop et al. 2003; Tonry 
et al. 2003; Riess et al. 2004). 

Unfortunately, the acoustic features in the matter cor- 
relations arc weak (10% contrast in the power spectrum) 
and on large scales. This means that one must survey very 
large volumes, of order lh~^ Gpc^, to detect the signature 
(Tegmark 1997; Goldberg & Strauss 1998; Eisenstein, Hu, 
& Tegmark 1998). Previous surveys have not found clean 
detections, due to sample size and (in some cases) survey 
geometry. Early surveys (Broadhurst et al. 1990; Landy et 
al. 1996; Einasto et al. 1997) found anomalous peaks that, 
though unlikely to be acoustic signatures (Eisenstein et al. 
1998), did not reappear in larger surveys. Percival et al. 

(2001) favored baryons at 2 cr in a power spectrum anal- 
ysis of data from the 2dF Galaxy Redshift Survey (here- 
after 2dFGRS; Colless et al. 2001), but Tegmark et al. 

(2002) did not recover the signal with 65% of the same 
data. Miller et al. (2002) argued that due to smearing 
from the window function, the 2dFGRS result could only 
be due to the excess of power on large scales in baryonic 
models and not due to a detection of the oscillations them- 
selves. A spherical harmonic analysis of the full 2dFGRS 
(Percival et al. 2004) did not find a significant baryon frac- 
tion. Miller et al. (2001) presented a "possible detection" 
(2.2 cr) from a combination of three smaller surveys. The 
analysis of the power spectrum of the SDSS Main sample 
(Tegmark et al. 2004a, b) did not address the question ex- 
plicitly but would not have been expected to detect the 
oscillations. The analysis of the correlation function of 
the 2dFGRS (Hawkins et al. 2003) and the SDSS Main 
sample (Zchavi ct al. 2004b) did not consider these scales. 
Large quasar surveys (Outram et al. 2003; Groom ct al. 
2004a, b; Yahata et al. 2004) are too limited by shot noise 
to reach the required precision, although the high redshift 
does give leverage on dark energy (Outram et al. 2004) via 
the Alcock-Paczynski (1979) test. 

In this paper, we present the large-scale correlation func- 
tion of a large spectroscopic sample of luminous, red galax- 
ies (LRGs) from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000). This sample covers 3816 square degrees out to 
a redshift of z = 0.47 with 46,748 galaxies. While it con- 
tains fewer galaxies than the 2dFGRS or the Main sample 
of the SDSS, the LRG sample (Eisenstein et al. 2001) has 
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been optimized for the study of structure on the largest 
scales and as a result it is expected to significantly out- 
perform those samples. First results on intermediate and 
small-scale clustering of the LRG sample were presented 
by Zchavi ct al. (2004a) and Eiscnstcin ct al. (2004), and 
the sample is also being used for the study of galaxy clus- 
ters and for the evolution of massive galaxies (Eisenstein 
et al. 2003). Here we focus on large scales and present the 
first clear detection of the acoustic peak at late times. 

The outline of the paper is as follows. We introduce 
the SDSS and the LRG sample in §2. In §3, we present 
the LRG correlation function and its covariance matrix, 
along with a discussion of tests we have performed. In §4, 
we fit the correlation function to theoretical models and 
construct measurements on the cosmological distance scale 
and cosmological parameters. We conclude in §5 with a 
general discussion of the results. Readers wishing to focus 
on the results rather than the details of the measurement 
may want to skip § 3.2, § 3.3, and § 4.2. 

2. THE SDSS LUMINOUS RED GALAXY SAMPLE 

The SDSS (York et al. 2000; Stoughton et al. 2002; 
Abazajian et al. 2003, 2004) is imaging 10^ square de- 
grees of high Galactic latitude sky in five passbands, u, g, 
r, i, and z (Fukugita et al. 1996; Gunn et al. 1998). Im- 
age processing (Lupton et al. 2001; Stoughton ct al. 2002; 
Pier et al. 2003; Ivezic et al. 2004) and calibration (Hogg 
et al. 2001; Smith, Tucker et al. 2002) allow one to se- 
lect galaxies, quasars, and stars for follow-up spectroscopy 
with twin fiber- fed doublc-spectrographs. Targets are as- 
signed to plug plates with a tiling algorithm that ensures 
nearly complete samples (Blanton ct al. 2003a): observing 
each plate generates 640 spectra covering 3800A to 9200A 
with a resolution of 1800. 

We select galaxies for spectroscopy by two algorithms. 
The primary sample (Strauss et al. 2002), referred to here 
as the SDSS Main sample, targets galaxies brighter than 
r = 17.77. The surface density of such galaxies is about 90 
per square degree, and the median redshift is 0.10 with a 
tail out to 2 ~ 0.25. The LRG algorithm (Eisenstein et al. 
2001) selects ~ 12 additional galaxies per square degree, 
using color-magnitude cuts in g, r, and i to select galaxies 
to a Petrosian magnitude r < 19.5 that are likely to be 
luminous early-types at redshifts up to ^ 0.5. All fluxes 
are corrected for extinction (Schlegel et al. 1998) before 
use. The selection is extremely efficient, and the redshift 
success rate is very high. A few additional galaxies (3 per 
square degree at z > 0.16) matching the rest-frame color 
and luminosity properties of the LRGs are extracted from 
the SDSS Main sample; we refer to this combined set as 
the LRG sample. 

For our clustering analysis, we use 46,748 luminous rod 
galaxies over 3816 square degrees and in the redshift range 
0.16 to 0.47. The sky coverage of the sample is shown in 
Hogg et al. (2004) and is similar to that of SDSS Data Re- 
lease 3 (Abazajian et al. 2004). Wo require that the galax- 
ies have rest-frame g-haiid absolute magnitudes —23.2 < 
Mg < -21.2 {h = 1, Ho = lOO/i km s'^ Mpc'^) where 
we have applied k corrections and passively evolved the 
galaxies to a fiducial redshift of 0.3. The resulting comov- 
ing number density is close to constant out to z = 0.36 
(i.e., volume limited) and drops thereafter; see Figure 1 



of Zehavi et al. (2004a). The LRG sample is unusual as 
compared to flux-limited surveys because it uses the same 
type of galaxy (luminous early- types) at all redshifts. 

We model the radial and angular selection functions us- 
ing the methods in the appendix of Zehavi ct al. (2004a). 
In brief, we build a model of the redshift distribution of 
the sample by integrating an empirical model of the lumi- 
nosity function and color distribution of the LRGs with 
respect to the luminosity-color selection boundaries of the 
sample. The model is smooth on small scales and includes 
the subtle interplay of the color-redshift relation and the 
color selection boundaries. To include slow evolutionary 
effects, we force the model to the observed redshift his- 
togram using low-pass filtering. We will show in § 3.2 that 
this filtering does not affect the correlation function on the 
scales of interest. 

The angular selection function is based on the spherical 
polygon description of lss_samplel4 (Blanton et al. 2004). 
We model fiber collisions and unobserved plates by the 
methods of Zehavi et al. (2004a). Regions with complete- 
ness below 60% (due to unobserved plates) are dropped 
entirely; the 3816 square degrees we use are highly com- 
plete. The survey mask excludes regions around bright 
stars, but does not otherwise model small-scale imperfec- 
tions in the data. 

The typical redshift of the sample is ^ = 0.35. The co- 
ordinate distance to this redshift is 962h~^ Mpc for flm = 
0.3, = 0.7. At this distance, lOOhr^ Mpc subtends 6.0° 
on the sky and corresponds to 0.04 in redshift. The LRG 
correlation length of 10h~^ Mpc (Zehavi et al. 2004a) sub- 
tends only 40' on the sky at z = 0.35; viewed on this scale, 
the survey is far more bulk than boundary. 

A common way to assess the statistical reach of a sur- 
vey, including the effects of shot noise, is by the effective 
volume (Feldman, Kaiser, & Peacock 1994; Tegmark 1997) 

where n{r) is the comoving number density of the sample 
at every location r. The effective volume is a hmction of 
the wavenumber k via the power amplitude P. For P = 
lQ*h-^Mpc^ {k « 0.15/1 Mpc"^), we find O.lS/j-^ Gpc''; 
for P = 4x IQ'^h-^ Mpc^ (fc w 0.05/?, Mpc"^), t)Mh-^ Gpc^; 
and for P = lO^/i-^Mpc^ (A; « 0.02/?, Mpc"^), .hbh'^ G-pc^ 
The actual survey volume is 0.72/i,~^ Gpc^, so there are 
roughly 700 cubes of 100/i~^ Mpc size in the survey. The 
relative sparseness of the LRG sample, n - IQ-^hl^ Mpc'^ 
is well suited to measuring power on large scales (Kaiser 
1986). 

In Figure 1, we show a comparison of the effective vol- 
ume of the SDSS LRG survey to other published surveys. 
While this calculation is necessarily a rough (~^30%) pre- 
dictor of statistical performance due to neglect of the exact 
survey boundary and our detailed assumptions about the 
amplitude of the power spectrum and the number density 
of objects for each survey, the SDSS LRG is clearly the 
largest survey to date for studying the linear regime by a 
factor of ~ 4. The LRG sample should therefore outper- 
form these surveys by a factor of 2 in fractional errors on 
large scales. Note that quasar surveys cover much more 
volume than even the LRG survey, but their effective vol- 
umes are worse, even on large scales, due to shot noise. 
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Fig. 1. — The effective volume (eq. [1]) as a function of wavenum- 
ber for various large redshift surveys. The effective volume is a rough 

— 1/2 

guide to the performance of a survey (errors scaling as V^^^ ) but 
should not be trusted to better than 30%. To facilitate comparison, 
we have assumed 3816 square degrees for the SDSS Main sample, 
the same area as the SDSS LRG sample presented in this paper and 
similar to the area in Data Release 3. This is about 50% larger than 
the sample analyzed in Tegmark et al. (2004a), which would be sim- 
ilar to the curve for the full 2dF Galaxy Redshift Survey (CoUess et 
al. 2003). We have neglected the potential gains on very large scales 
from the 99 outrigger fields of the 2dFGRS. The other surveys are 
the MX survey of clusters (Miller & Batuski 2001), the PSCz survey 
of galaxies (Sutherland et al. 1999), and the 2QZ survey of quasars 
(Groom et al. 2004a). The SDSS DR3 quasar survey (Schneider et 
al. 2005) is similar in effective volume to the 2QZ. For the amplitude 
of P{k), we have used ag = 1 for 2QZ and PSGz and 3.6 for the 
MX survey. We used erg = 1.8 for SDSS LRG, SDSS Main, and the 
2dFGRS; For the latter two, this value represents the amplitude of 
clustering of the luminous galaxies at the surveys' edge; at lower 
redshift, the number density is so high that the choice of trg is ir- 
relevant. Reducing SDSS Main or 2dFGRS to trg = 1, the value 
typical of normal galaxies, decreases their Veff by 30%. 



3. THE REDSHIFT-SPACE CORRELATION FUNCTION 

3.1. Correlation function estimation 

In this paper, we analyze the large-scale clustering us- 
ing the two-point correlation function (Peebles 1980). In 
recent years, the power spectrum has become the common 
choice on large scales, as the power in different Fourier 
modes of the linear density field is statistically indepen- 
dent in standard cosmology theories (Bardeen et al. 1986). 
However, this advantage breaks down on small scales due 
to non-linear structure formation, while on large scales, 
elaborate methods are required to recover the statistical 
independence in the face of survey boundary effects (for 
discussion, see Tegmark et al. 1998). The power spec- 
trum and correlation function contain the same informa- 
tion in principle, as they are Fourier transforms of one 
another. The property of the independence of different 
Fourier modes is not lost in real space, but rather it is en- 
coded into the off-diagonal elements of the covariance ma- 
trix via a linear basis transformation. One must therefore 
accurately track the full covariance matrix to use the cor- 
relation function properly, but this is feasible. An advan- 
tage of the correlation function is that, unlike in the power 
spectrum, small-scale effects like shot noise and intra-halo 
astrophysics stay on small scales, well separated from the 



linear regime fluctuations and acoustic effects. 

We compute the redshift-space correlation function us- 
ing the Landy-Szalay estimator (Landy & Szalay 1993). 
Random catalogs containing at least 16 times as many 
galaxies as the LRG sample were constructed according to 
the radial and angular selection functions described above. 
We assume a flat cosmology with il„i = 0.3 and SIa — 0.7 
when computing the correlation function. We place each 
data point in its comoving coordinate location based on 
its redshift and compute the comoving separation between 
two points using the vector difference. We use bins in sep- 
aration of 4/i~^Mpc from 10 to 30ft.~^Mpc and bins of 
lOh^^ Mpc thereafter out to 180/1"^ Mpc, for a total of 20 
bins. 

We weight the sample using a scale-independent weight- 
ing that depends on redshift. When computing the corre- 
lation function, each galaxy and random point is weighted 
by 1/(1 -t- n{z)Pyj) (Feldman, Kaiser, & Peacock 1994) 
where n{z) is the comoving number density and Pw = 
40, OOO/i^'^ Mpc^. We do not allow to change with scale 
so as to avoid scale-dependent changes in the effective bias 
caused by differential changes in the sample redshift. Our 
choice of P^ is close to optimal at fc w 0.05/iMpc~^ and 
within 5% of the optimal errors for all scales relevant to 
the acoustic oscillations (fc < 0.15/i Mpc~^). At z < 0.36, 
nPi^ is about 4, while nP^ « 1 at z = 0.47. Our results 
do not qualitatively depend upon the value of Pw ■ 

Redshift distortions cause the redshift-space correlation 
function to vary according to the angle between the sep- 
aration vector and the line of sight. To ease comparison 
to theory, we focus on the spherically averaged correla- 
tion function. Because of the boundary of the survey, the 
number of possible tangential separations is somewhat un- 
derrepresented compared to the number of possible line of 
sight separations, particularly at very large scales. To cor- 
rect for this, we compute the correlation functions in four 
angular bins. The effects of redshift distortions are ob- 
vious: large-separation correlations are smaller along the 
line of sight direction than along the tangential direction. 
We sum these four correlation functions in the proportions 
corresponding to the fraction of the sphere included in the 
angular bin, thereby recovering the spherically averaged 
redshift-space correlation function. We have not yet ex- 
plored the cosmological implications of the anisotropy of 
the correlation function (Matsubara & Szalay 2003). 

The resulting redshift-space correlation function is shown 
in Figure 2. A more convenient view is in Figure 3, where 
we have multiplied by the square of the separation, so as 
to flatten out the result. The errors and overlaid models 
will be discussed below. The bump at 100/i^^Mpc is the 
acoustic peak, to be described in §4.1. 

The clustering bias of LRGs is known to be a strong 
function of luminosity (Hogg et al. 2003; Eisenstein et al. 
2004; Zehavi et al. 2004a) and while the LRG sample is 
nearly volume-limited out to z ~ 0.36, the flux cut does 
produce a varying luminosity cut at higher redshifts. If 
larger scale correlations were preferentially drawn from 
higher redshift, we would have a differential bias (see dis- 
cussion in Tegmark et al. 2004a). However, Zehavi et al. 
(2004a) have studied the clustering amplitude in the two 
limiting cases, namely the luminosity threshold at z < 0.36 
and that at z = 0.47. The differential bias between these 
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Fig. 2. — The large-scale redsliift-space correlation function of the 
SDSS LRG sample. The error bars are from the diagonal elements 
of the mock-catalog covariance matrix; however, the points are cor- 
related. Note that the vertical axis mixes logarithmic and linear 
scalings. The inset shows an expanded view with a linear vertical 
axis. The models are Umh'^ = 0.12 (top, green), 0.13 (red), and 
0.14 (bottom with peak, blue), all with Q^/i^ = 0.024 and n = 0.98 
and with a mild non-linear prescription folded in. The magenta 
line shows a pure CDM model {Qmh^ = 0.105), which lacks the 
acoustic peak. It is interesting to note that although the data ap- 
pears higher than the models, the covariance between the points is 
soft as regards overall shifts in ^(s). Subtracting 0.002 from ^{s) 
at all scales makes the plot look cosmetically perfect, but changes 
the best-fit by only 1.3. The bump at 100/i~^ Mpc scale, on the 
other hand, is statistically significant. 

two samples on large scales is modest, only 15%. We make 
a simple parameterization of the bias as a function of red- 
shift and then compute averaged as a function of scale 
over the pair counts in the random catalog. The bias varies 
by less than 0.5% as a function of scale, and so we conclude 
that there is no effect of a possible correlation of scale with 
redshift. This test also shows that the mean redshift as a 
function of scale changes so little that variations in the 
clustering amplitude at fixed luminosity as a function of 
redshift are negligible. 

3.2. Tests for systematic errors 

We have performed a number of tests searching for po- 
tential systematic errors in our correlation function. First, 
we have tested that the radial selection function is not in- 
troducing features into the correlation function. Our selec- 
tion function involves smoothing the observed histogram 
with a box-car smoothing of width Az = 0.07. This cor- 
responds to reducing power in the purely radial mode at 
k = 0.03/1 Mpc"^ by 50%. Purely radial power at k = 0.04 
(0.02)/iMpc"i is reduced by 13% (86%). The effect of this 
suppression is negligible, only 5 x lO^'* (10"'*) on the cor- 
relation function at the 30 (100) h^^Mpc scale. Simply 
put, purely radial modes are a small fraction of the total 
at these wavelengths. We find that an alternative radial 
selection function, in which the redshifts of the random 
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Fig. 3. — As Figure 2, but plotting the correlation function times 
s^. This shows the variation of the peak at 20h~^ Mpc scales that is 
controlled by the redshift of equality (and hence by Q,mh?). Vary- 
ing ^Imh'^ alters the amount of large-to-small scale correlation, but 
boosting the large-scale correlations too much causes an inconsis- 
tency at 30h~^ Mpc. The pure CDM model (magenta) is actually 
close to the best-fit due to the data points on intermediate scales. 



catalog are simply picked randomly from the observed red- 
shifts, produces a negligible change in the correlation func- 
tion. This of course corresponds to complete suppression 
of purely radial modes. 

The selection of LRGs is highly sensitive to errors in the 
photometric calibration of the g, r, and i bands (Eisenstein 
et al. 2001). We assess these by making a detailed model 
of the distribution in color and luminosity of the sample, 
including photometric errors, and then computing the vari- 
ation of the number of galaxies accepted at each redshift 
with small variations in the LRG sample cuts. A 1% shift 
in the r — i color makes a 8-10% change in number den- 
sity; a 1% shift in the g — r color makes a 5% changes in 
number density out to z = 0.41, dropping thereafter; and 
a 1% change in all magnitudes together changes the num- 
ber density by 2% out to z = 0.36, increasing to 3.6% at 
z = 0.47. These variations are consistent with the changes 
in the observed redshift distribution when we move the 
selection boundaries to restrict the sample. Such photo- 
metric calibration errors would cause anomalies in the cor- 
relation function as the square of the number density vari- 
ations, as this noise source is uncorrelated with the true 
sky distribution of LRGs. 

Assessments of calibration errors based on the color of 
the stellar locus find only 1% scatter in g, r, and i (Ivezic 
et al. 2004), which would translate to about 0.02 in the 
correlation function. However, the situation is more favor- 
able, because the coherence scale of the calibration errors 
is limited by the fact that the SDSS is calibrated in regions 
about 0.6° wide and up to 15° long. This means that there 
are 20 independent calibrations being applied to a given 
6° (100ft.~^Mpc) radius circular region. Moreover, some 
of the calibration errors are even more localized, being 
caused by small mischaracterizations of the point spread 
function and errors in the flat field vectors early in the 
survey (Stoughton et al. 2002). Such errors will average 
down on larger scales even more quickly. 

The photometric calibration of the SDSS has evolved 
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slightly over time (Abazajian et al. 2003), but of course the 
LRG selection was based on the calibrations at the time 
of targeting. We make our absolute magnitude cut on the 
latest ('uber') calibrations (Blanton et al. 2004; Finkbciner 
et al. 2004), although this is only important at z < 0.36 
as the targeting flux cut limits the sample at higher red- 
shift. We test whether changes in the calibrations might 
alter the correlation function by creating a new random 
catalog, in which the differences in each band between the 
target epoch photometry and the 'uber'-calibration values 
are mapped to angularly dependent redshift distributions 
using the number density derivatives presented above. Us- 
ing this random catalog makes negligible differences in the 
correlation function on all scales. This is not surprising: 
while an rms error of 2% in r — i would give rise to a 
0.04 excess in the correlation function, scales large enough 
that this would matter have many independent calibra- 
tions contributing to them. On the other hand, the 'uber'- 
calibration of the survey does not necessarily fix all cali- 
bration problems, particularly variations within a single 
survey run, and so this test cannot rule out an arbitrary 
calibration problem. 

We continue our search for calibration errors by breaking 
the survey into 10 radial pieces and measuring the cross- 
correlations between the non-adjacent slabs. Calibration 
errors would produce significant correlations on large an- 
gular scales. Some cross-correlations have amplitudes of 
2-3%, but many others don't, suggesting that this is sim- 
ply noise. We also take the full matrix of cross-correlations 
at a given separation and attempt to model it (minus the 
diagonal and first off-diagonal elements) as an outer prod- 
uct of vector with itself, as would be appropriate if it were 
dominated by a single type of radial perturbation, but we 
do not find plausible or stable vectors, again indicative of 
noise. Hence, we conclude that systematic errors in ^(r) 
due to calibration must be below 0.01. 

It is important to note that calibration errors in the 
SDSS produce large-angle correlations only along the scan 
direction. Even if errors were noticeably large, they would 
not produce narrow features such as that seen at the 100/i~^ 
Mpc scale for two reasons. First, the projections from the 
three-dimensional sphere to one-dimension strips on the 
sky necessarily means that a given angular scale maps to 
a wide range of three-dimensional separations. Second, the 
comoving angular diameter distance used to translate an- 
gles into transverse separations varies by a factor of three 
from z = 0.16 to 0.47, so that a preferred angle would not 
map to a narrow range of physical scales. We therefore 
expect that calibration errors would appear as a smooth 
anomalous correlation, rolling off towards large scales. 

Breaking the sample into two redshift slices above and 
below z — 0.36 yields similar correlation functions (Fig. 4). 
Errors in calibration or in the radial selection function 
would likely enter the two redshift slices in different man- 
ners, but we see no sign of this. In particular, the bump 
in the correlation function appears in both slices. While 
this could in principle give additional leverage on the cos- 
mological distance scale, we have not pursued this in this 
paper. 
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Fig. 4. — The correlation function for two different redshift shoes, 
0.16 < z < 0.36 (filled squares, black) and 0.36 < z < 0.47 (open 
squares, red). The latter is somewhat noisier, but the two are quite 
similar and both show evidence for the acoustic peak. Note that the 
vertical axis mixes logarithmic and linear scalings. 



3.3. Covariance Matrix 

Because of the large number of separation bins and the 
large scales being studied, it is infeasible to use jackknife 
sampling to construct a covariance matrix for ^(r). In- 
stead, we use a large set of mock catalogs to construct a 
covariance matrix and then test that matrix with a smaller 
number of jackknife samples. 

Our mock catalogs are constructed using PTHalos (Scoc- 
cimarro & Sheth 2002) with a halo occupation model (Ma 
& Fry 2000; Seljak 2000; Peacock & Smith 2000; Scocci- 
marro et al. 2001; Cooray & Sheth 2002; Berlind et al. 
2003) that matches the observed amplitude of clustering 
of LRGs (Zchavi et al. 2004a). PTHalos distributes dark 
matter halos according to extended Press-Schechter theory 
conditioned by the large-scale density field from second- 
order Lagrangian perturbation theory. This generates den- 
sity distributions that fully include Gaussian linear the- 
ory, but also include second-order clustering, redshift dis- 
tortions, and the small-scale halo structure that domi- 
nates the non-Gaussian signal. We use a cosmology of 
fim = 0.3, JIa = 0.7, h = 0.7 for the mock catalogs. We 
subsample the catalogs so as to match the comoving den- 
sity of LRGs as a function of redshift exactly, but do not 
attempt to include the small redshift dependence in the 
amplitude of the bias. Our catalogs match the angular 
geometry of the survey except in some fine details involv- 
ing less than 1% of the area. Two simulation boxes, each 
(1250/1-1 Mpc)2 X 2500/1-1 Mpc, one for the North Galac- 
tic Cap and the other for the South, are pasted together in 
each catalog. These two regions of the survey are very well 
separated in space, so this division is harmless. We gen- 
erate 1278 mock catalogs with independent initial condi- 
tions, compute the correlation function in each, and com- 
pute the covariance matrix from the variations between 
them. 

The resulting matrix shows considerable correlations be- 
tween neighboring bins, but with an enhanced diagonal 
due to shot noise. The power on scales below 10h~^ Mpc 
creates significant correlation between neighboring scales. 
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A curious aspect is that on large scales, where our bins are 
10h~^ Mpc wide, the reduced inverse covariance matrix is 
quite close to tri-diagonal; the first off-diagonal is typically 
around 0.4, but the second and subsequent are typically 
a few percent. Such matrices correspond to exponential 
decays of the off-diagonal correlations (Rybicki & Press 
1994). 

We test the covariance matrix in several ways. First, 
as described in §4.3, the best-fit cosmological model has 
= 16.1 on 17 degrees of freedom {p = 0.52), indicat- 
ing that the covariances are of the correct scale. Next, we 
subdivide the survey into ten large, compact subregions 
and look at the variations between jack-knifed samples 
(i.e., computing ^ while excluding one region at a time). 
The variations of ^(r) between the ten jackknife samples 
matches the diagonal elements of the covariance matrix 
to within 10%. We then test the off-diagonal terms by 
asking whether the jack-knife residuals vary appropriately 
relative to the covariance matrix. If one defines the resid- 
ual A^j between the jth jackknife sample and the mean, 

then the construction "J^ab ^0(^a)'-^ai)^^^j where C is 
the covariance matrix and the sums are over the 20 radial 
bins, should have a mean (averaged over the ten jackknife 
sample) of about 20/(10 — 1) « 2.2. The mean value is 
2.7 ± 0.5, in reasonaljle agreement. 

We get similar results for cosmological parameters when 
using a covariance matrix that is based on constructing 
the Gaussian approximation (Feldman, Kaiser, & Pea- 
cock 1994; Tegmark 1997) of independent modes in Fourier 
space. We use the effective volume at each wavenumber, 
plus an extra shot noise term to represent non-Gaussian 
halos, and then rotate this matrix from the diagonal Fourier 
basis into the real-space basis. This matrix gives reason- 
able values, satisfies the jackknife tests, and gives sim- 
ilar values and errors on cosmological parameters. 

Finally, wc test our results for cosmological parameters 
by using the following hybrid scheme. We use the mock- 
catalog covariance matrix to find the best-fit cosmological 
model for each of the ten jackknife samples, and then use 
the rms of the ten best-fit parameter sets to determine 
the errors. This means that we are using the mock cata- 
logs to weight the ^(r) measurements, but relying on the 
jackknifed samples to actually determine the variance on 
the cosmological parameters. We will qiiote the numerical 
results in § 4.3, but here we note that the resulting con- 
straints on the acoustic scale and matter density match 
the errors inferred from the fitting with the mock-catalog 
covariance matrix. Hence, we conclude that our covariance 
matrix is generating correct results for our model fitting 
and cosmological parameter estimates. 

4. CONSTRAINTS ON COSMOLOGICAL MODELS 

4.1. Linear theory 

Given the value of the matter density ilmh"^, the baryon 
density flbh^, the spectral tilt n, and a possible neutrino 
mass, adiabatic cold dark matter (CDM) models predict 
the linear matter power spectrum (and correlation func- 
tion) up to an amplitude factor. There are two primary 
physical scales at work (for discussions, see Hu et al. 1997; 
Eisenstein & Hu 1998; Hu & Dodelson 2002). First, the 
clustering of CDM is suppressed on scales small enough 
to have been traversed by the neutrinos and/or photons 



during the radiation-dominated period of the universe. 
This introduces the characteristic turnover in the CDM 
power spectrum at the scale of the horizon at matter- 
radiation equality. This length scales as {^l„ih'^)~^ (assum- 
ing the standard cosmic neutrino background). Second, 
the acoustic oscillations have a characteristic scale known 
as the sound horizon, which is the comoving distance that 
a sound wave can travel between the big bang and recom- 
bination. This depends both on the expansion history of 
the early universe and on the sound speed in the plasma; 
for models close to the concordance cosmology, this length 
scales as (Q„/i2)-0-25(O(,/i2)-0.08 ^jju 2004). The spectral 
tilt and massive neutrinos tilt the low-redshift power spec- 
trum (Bond & Szalay 1983), but don't otherwise change 
these two scales. Importantly, given ftmh^ and flbh^, the 
two physical scales are known in absolute units, with no 
factors of H^^. 

Wc use CMBfast (Seljak & Zaldarriaga 1996; Zaldar- 
riaga et al. 1998; Zaldarriaga & Seljak 2000) to compute 
the linear power spectra, which we convert to correlation 
functions with a Fourier transform. It is important to note 
that the harmonic series of acoustic peaks found in the 
power spectrum transform to a single peak in the correla- 
tion function (Matsubara 2004). The decreasing envelope 
of the higher harmonics, due to Silk damping (Silk 1968) 
or non-linear gravity (Meiksin et al. 1999), corresponds to 
a broadening of the single peak. This broadening decreases 
the accuracy with which one can measure the ccntroid of 
the peak, equivalent to the degradation of acoustic scale 
measurements caused by the disappearance of higher har- 
monics in the power spectrum (Sco & Eisenstein 2003). 

Examples of the model correlation functions arc shown 
in Figures 2 and 3, each with fib/i^ = 0.024 and h = 0.7 
but with three values of ^mh^- Higher values of ^mh^ 
correspond to earlier epochs of matter-radiation equality, 
which increase the amount of small-scale power compared 
to large, which in turn decreases the correlations on large 
scales when holding the small-scale amplitude fixed. The 
acoustic scale, on the other hand, depends only weakly on 

4.2. Non-linear corrections 

The precision of the LRG correlations is such that we 
cannot rely entirely on linear theory even at r > lOft."^ Mpc 
(wavenumbers fc < O.l/iMpc"^). Non-linear gravity (Meiksin 
et al. 1999), redshift distortions (Kaiser 1987; Hamilton 
1998; Scoccimarro 2004), and scale-dependent bias all en- 
ter at a subtle level. We address these by applying correc- 
tions derived from N-body simulations and the Smith et 
al. (2003) halof it method. 

Seo & Eisenstein (2005) present a set of 51 N-body sim- 
ulations of the WMAP best-fit cosmology (Spergel et al. 
2003). Each simulation has 256^ particles and is 512h~^ Mpc 
in size. The outputs at 2; = 0.3 were analyzed to find their 
power spectra and correlation functions, in both real and 
redshift-space, including simple halo bias. This provides 
an accurate description of the non- linear gravitational and 
redshift distortion effects on large scales. 

For a general cosmology, we begin from the CMBfast 
linear power spectrum. We next correct for the erasure of 
the higher acoustic peak that occurs due to mode-coupling 
(Meiksin et al. 1999). We do this by generating the "no- 
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Fig. 5. — Scale-dependent corrections derived from 51 N-body 
simulations, each 512h~^ Mpc comoving with 256"^ particles (Seo & 
Eisenstein 2005). The crosses show the ratio between the non-linear 
matter correlation function and the linear correlation function; the 
dashed line is the model we use from Smith et al. (2003). The 
solid points are the ratio between the biased correlation function 
(using a simple halo mass cut) to the non-linear matter correlation 
function. The open squares are the ratio of the biased redshift-space 
correlation function to the biased real-space correlation function, 
after removing the large-scale asymptotic value (Kaiser 1987), which 
we simply fold into the correlation amplitude parameter. The open 
triangles show the product of these two effects, and the solid line 
is our fit to this product. These corrections are of order 10% at 
lOh^^ Mpc separations and decrease quickly on larger scales. In 
addition to these corrections, we mimic the erasure of the small-scale 
acoustic oscillations in the power spectrum by using a smoothed 
cross-over at fe = 0.14/iMpc~^ between the CMBfast linear power 
spectrum and the no-wiggle form from Eisenstein & Hu (1998). 



wiggle" approximation from Eisenstein & Hu (1998), which 
matches the overall shape of the linear power spectrum 
but with the acoustic oscillations edited out, and then 
smoothly interpolating between the linear spectrum and 
the approximate one. We use a Gaussian exp[— (fca)^], 
with a scale a = 7h^^ Mpc chosen to approximately match 
the suppression of the oscillations seen in the power spec- 
trum of the N-body simulations. 

We next use the Smith et al. (2003) package to compute 
the alterations in power from non-linear gravitational col- 
lapse. We use a model with F = 0.162 and ag = 0.85. 
This model has the same shape for the power spectrum 
on small scales as the WMAP best-fit cosmology and is 
therefore a reasonable zcro-baryon model with which to 
compute. We find the quotient of the non-linear and lin- 
ear power spectra for this model and multiply that ratio 
onto our baryonic power spectrum. We then Fourier trans- 
form the power spectrum to generate the real-space corre- 
lation function. The accuracy of this correction is shown 
in the bottom of Figure 5: the crosses are the ratio in the 
N-body simulations between the z — 0.3 and z — 49 corre- 
lation functions, while the dashed line is the Smith et al. 
(2003) derived correction. 

Next we correct for redshift distortions and halo bias. 
We use the N-body simulations to find the ratio of the 
redshift-space biased correlation function to the real-space 
matter correlation function for a halo mass threshold that 
approximately matches the observed LRG clustering am- 
plitude. This ratio approaches an asymptotic value on 



large scales that we simply include into the large-scale 
clustering bias. After removing this asymptotic value, the 
remaining piece of the ratio is only 10% at r = 10h~^ Mpc. 
We fit this to a simple smooth function and then multiply 
the model correlation functions by this fit. Figure 5 shows 
the N-body results for the ratio of redshift-space ^ to real- 
space ^ for the halos and the ratio of the real-space halo 
correlation function to that of the matter. The solid line 
shows our fit to the product. Tripling the mass thresh- 
old of the bias model increases the correction by only 30% 
(i.e., 1.12 to 1.16 on small scales). 

All of the corrections in Figure 5 are clearly small, only 
10% at r > lOh-^ Mpc. This is because lO/i"^ Mpc sep- 
arations are larger than any virialized halo and arc only 
affected by the extremes of the finger of God redshift dis- 
tortions. While our methods are not perfect, they are 
plausibly matched to the allowed cosmology and to the 
bias of the LRGs. As such, we believe that the corrections 
should be accurate to a few percent, which is sufficient for 
our purposes. Similarly, while we have derived the correc- 
tions for a single cosmological model, the data constrain 
the allowed cosmologies enough that variations in the cor- 
rections will be smaller than our tolerances. For example, 
increasing as to 1.0 for the halofit calculation changes the 
corrections at r > 10h~^ Mpc by less than 2%. 

We stress that while galaxy clustering bias does rou- 
tinely affect large-scale clustering (obviously so in the LRG 
sample, with bias w 2), it is very implausible that it 
would mimic the acoustic signature, as this would require 
galaxy formation physics to have a strong preferred scale 
at I00h~^ Mpc. Galaxy formation prescriptions that in- 
volve only small-scale physics, such as that involving dark 
matter halos or even moderate scale radiation transport, 
necessarily produce smooth effects on large scales (Coles 
1993; Pry & Gaztanaga 1993; Scherrer & Weinberg 1998). 
Even long-range effects that might be invoked would need 
to affect lOOft."^ Mpc scales differently from 80 or 130. Our 
detection of the acoustic peak cannot reasonably be ex- 
plained as an illusion of galaxy formation physics. 

4.3. Measurements of the acoustic and equality scales 

The observed LRG correlation function could differ from 
that of the correct cosmological model in amplitude, be- 
cause of clustering bias and uncertain growth functions, 
and in scale, because we may have used an incorrect cos- 
mology in converting from redshift into distance. Our goal 
is to use the comparison between observations and theory 
to infer the correct distance scale. 

Note that in principle a change in the cosmological model 
would change the distances differently for different red- 
shifts, requiring us to recompute the correlation function 
for each model choice. In practice, the changes are small 
enough and the redshifts close enough that we treat the 
variation as a single dilation in scale (similar to Blake & 
Glazebrook 2003). This would be a superb approximation 
at low redshift, where all distances behave inversely with 
the Hubble constant. By z = 0.35, the effects of cosmolog- 
ical acceleration are beginning to enter. However, we have 
checked explicitly that our single-scale approximation is 
good enough for between 0.2 and 0.4. Relative to our 
fiducial scale at z = 0.35, the change in distance across the 
redshift range 0.16 < z < 0.47 is only 3% peak to peak 
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Fig. 6. — The values of the models as a function of the dilation 
of the scale of the correlation function. This corresponds to altering 
Dv{0.35) relative to the baseline cosmology of f2 = 0.3, Qa = 0.7, 
h = 0.7. Each line (save the magenta line) in the plot is a different 
value of Qmh^, 0.11, 0.13, and 0.15 from left to right. Qth^ = 0.024 
and n = 0.98 are used in all cases. The amplitude of the model 
has been marginalized over. The best-fit is 16.1 on 17 de grees 
of freedom, consistent with expectations. The magenta line (open 
symbols) shows the pure CDM model with Qmh'^ = 0.10; it has a 
best of 27.8, which is rejected at 3.4 a. Note that this curve 
is also much broader, indicating that the lack of an acoustic peak 
makes the scale less constrainable. 



for flm = 0.2 compared to 0.3, and even these variations 
largely cancel around the z = 0.35 midpoint where we will 
quote our cosmological constraints. 

The other error in our one-scale-parameter approxima- 
tion is to treat the line-of-sight dilation equivalently to 
the transverse dilation. In truth, the Hubble parame- 
ter changes differently from the angular diameter distance 
(the Alcock-Paczynski (1979) effect). For small deviations 
from Q„i ~ 0.3, flA — 0.7, the change in the Hubble pa- 
rameter at z — 0.35 is about half of that of the angular 
diameter distance. We model this by treating the dilation 
scale as the cube root of the product of the radial dilation 
times the square of the transverse dilation. In other words, 
we define 

ll/3 



Dviz) 



Dm{zY 



H{z) 



(2) 



where II{z) is the Hubble parameter and Dm{z) is the 
comoving angular diameter distance. As the typical red- 
shift of the sample is z = 0.35, we quote our result for the 
dilation scale as I?y(0.35). For our fiducial cosmology of 
= 0.3, f^A = 0.7, h = 0.7, DviO.ib) = 1334 Mpc. 
We compute parameter constraints by computing 
(using the full covariance matrix) for a grid of cosmological 
models. In addition to cosmological parameters of flmh^, 
flbh^, and n, we include the distance scale Z?y (0.35) of the 
LRG sample and marginalize over the amplitude of the 
correlation function. Parameters such as h, Qm, ^k, and 
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Fig. 7. — The likelihood contours of CDM models as a func- 
tion of Qmh^ and Z)\/(0.35). The likelihood has been taken to be 
proportional to exp(— x^/2), and contours corresponding to 1 c 
through 5 a for a 2-d Gaussian have been plotted. The one- 
dimensional marginalized values are Umh'^ = 0.130 ± 0.010 and 
Dv(0.35) = 1370 ± 64 Mpc. We overplot lines depicting the two 
major degeneracy directions. The solid (red) line is a line of con- 
stant Clmh^DY{0.35), which would be the degeneracy direction for 
a pure CDM model. The dashed (magenta) line is a line of con- 
stant sound horizon, holding Q^h^ = 0.024. The contours clearly 
deviate from the pure CDM degeneracy, implying that the peak at 
100/i~^ Mpc is constraining the fits. 



w{z) are subsumed within Dy(0.35). We assume h = 0.7 
when computing the scale at which to apply the non-linear 
corrections; having set those corrections, we then dilate the 
scale of the final correlation function. 

The WMAP data (Bennett et al. 2003), as well as com- 
binations of WMAP with large-scale structure (Spergel et 
al. 2003; Tegmark et al. 2004b), the Lyman-alpha forest 
(McDonald et al. 2004; Seljak et al. 2004), and big bang 
nucleosynthesis (e.g., Buries et al. 2001; Coc et al. 2004), 
constrain fl(,h^ and n rather well and so to begin, we hold 
these parameters fixed (at 0.024 and 0.98, respectively), 
and consider only variations in Vlmh?- In practice, the 
sound horizon varies only as (ri^/i^)^'^'^^, which means 
that the tight constraints from WMAP (Spergel et al. 
2003) and big bang nucleosynthesis (Buries et al. 2001) 
make the uncertainties in the baryon density negligible. 

Figure 6 shows as a function of the dilation for three 
different values of Q.mh^ , 0.11, 0.13, and 0.15. Scanning 
across all ilmh^, the best-fit is 16.1 on 17 degrees of free- 
dom (20 data points and 3 parameters: flmh'^, Z)y(0.35), 
and the amplitude). Figure 7 shows the contours of equal 
in flmh^ and £'y(0.35), corresponding to 1 cr up to 5 cr 
for a 2-dimensional Gaussian likelihood function. Adopt- 
ing a hkelihood proportional to exp(— x^/2), we project 
the axes to find fl^h'^ = 0.130 ± 0.010 and 15^(0.35) = 
1370 ± 64Mpc (4.7%), where these are 1 cr errors. 

Figure 7 also contains two lines that depict the two phys- 
ical scales. The solid line is that of constant ^Imh^Dy, 
which would place the (matter-radiation) equality scale at 
a constant apparent location. This would be the degener- 
acy direction for a pure cold dark matter cosmology and 
would be a line of constant F = Qmh were the LRG sam- 
ple at lower redshift. The dashed line holds constant the 
sound horizon divided by the distance, which is the appar- 
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Table 1 



Summary of Parameter Constraints from LRGs 

1.2 



1700 



Dv(0.35) 

Ro.35 = Dv(0.3 5)/Dai (1089) 
A = D^^(0.35)y''n„H2/0.35c 0.469(n/0.98)~"-35 _|_ 0.017 (3.6%) 



0.130(n/0.98)i-2 ±0.011 
1370 ± 64 Mpc (4.7%) 
0.0979 ± 0.0036 (3.7%) 



NOTES. — We assume Qi,h^ = 0.024 throughout, but variations per- 
mitted by WMAP create negligible changes here. We use n = 0.98, 
but where variations by 0.1 would create 1 cr changes, we include an 
approximate dependence. The quantity A is discussed in § 4.5. All 
constraints are 1 a. 



ent location of the acoustic scale. One sees that the long 
axis of the contours falls in between these two, and the fact 
that neither direction is degenerate means that both the 
equality scale and the acoustic scale have been detected. 
Note that no information from the CMB on flmh^ has been 
used in computing x^, and so our constraint on Q.m.h'^ is 
separate from that from the CMB. 

The best-fit pure CDM model has = 27.8, which 
means that it is disfavored by Ax^ = 11-7 compared to 
the model with Q.^h'^ = 0.024. Note that we are not 
marginalizing over the baryon density, so these two pa- 
rameter spaces have the same number of parameters. The 
baryon signature is therefore detected at 3.4 a. As a more 
stringent version of this, we find that the baryon model 
is preferred by Ax^ — 8.8 (3.0 cr) even if we only include 
data points between 60 and 180ft.~^Mpc. Figure 6 also 
shows x^ for the pure-CDM model as a function of di- 
lation scale; one sees that the scale constraint on such a 
model is a factor of two worse than the baryonic models. 
This demonstrates the importance of the acoustic scale in 
our distance inferences. 

As most of our distance leverage is coming from the 
acoustic scale, the most robust distance measurement we 
can quote is the ratio of the distance to z = 0.35 to the 
distance to z = 1089 (the redshift of decoupling, Bennett 
et al. 2003). This marginalizes over the uncertainties in 
flmh^ and would cancel out more exotic errors in the sound 
horizon, such as from extra relativistic species (Eisenstein 
& White 2004). We denote this ratio as 

_ i?v(0.35) 

Note that the CMB measures a purely transverse distance, 
while the LRG sample measures the hybrid in Equation 
(2). i:>A/(1089) = 13700 Mpc for the fl^ = 0.3, flA = 0.7, 
h = 0.7 cosmology (Tegmark et al. 2004b), with uncer- 
tainties due to imperfect measurement of the CMB an- 
gular acoustic scale being negligible at < 1%. We find 
-R0.35 = 0.0979 ± 0.0036, which is a 4% measurement of 
the relative distance to z = 0.35 and z — 1089. Table 1 
summarizes our numerical results on these basic measure- 
ments. 

To stress that the acoustic scale is responsible for the 
distance constraint, we repeat our fitting having discarded 
the two smallest separation bins (10 < s < 18h~^ Mpc) 
from the correlation function. This is shown in Figure 8. 
One sees that the constraints on ilmh^ have degraded (to 




0.14 



Fig. 8. — As Figure 7, but now with scales below 18h~^ Mpc ex- 
cluded from the computation. This leaves 18 separation bins and 
15 degrees of freedom. The contours are now obviously aligned to 
the line of constant sound horizon, and the constraints in the Qmh'^ 
direction are weakened by 40%. As Figure 3 would suggest, the data 
at scales below 18/i~^ Mpc help to constrain Q,mh? , twisting the 
contours towards the pure CDM degeneracy. Dropping the smaller 
scales doesn't affect the constraint on -R0.35; we find 0.0973 ± 0.0038 
as compared to 0.0979 ± 0.0036 before. 
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Fig. 9. — As Figure 7, but now with a spectral tilt of n = 0.90. The 
best fit has = 17.8. The primary effect is a shift to larger Qrnh'^, 
0.143 ± 0.011. However, this shift occurs at essentially constant 
we find 0.0986 ± 0.0041. Again, the acoustic scale robustly 
determines the distance, even though the spectral tilt biases the 
measurement of the equality scale. 



0.136±0.014), but the contours remain well confined along 
the direction of constant acoustic scale (the dashed line) . 
We find a distance ratio of 0.0973 ± 0.0038, essentially 
identical to what we found above, with a best x^ of 13.7 
on 15 degrees of freedom. 

Varying the spectral tilt, which has a similar effect to 
including massive neutrinos, is partially degenerate with 
n„^h^, but the ratio of the distances is very stable across 
the plausible range. Repeating the fitting with n — 0.90 
changes the distance ratio i?o.35 to 0.0986 ± 0.0041, a less 
than 1% change. The distance itself changes by only 2% to 
Dv{0.35) = 1344 ± 70. The change in the likehhood con- 
tours is shown in Figure 9. The best-fit Qmh"^ for n = 0.90 
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Fig. 10. — a) As Figure 7, but overplotted with model predictions from constant w flat models. For a given value of Qmh^ and w, the 
angular scale of the CMB acoustic peaks (known to 1%) determines Clm and Hq. Of course, the required fim is a function of w and Q,mh? ■ 
The solid red lines show lines of constant w; the dashed lines show lines of constant Q.m- Our knowledge of Q,mh? still limits our inference of 
w. b) As (a), but the dashed lines are now lines of constant Hq. 



is 0.143 ± 0.011, and so we approximate our ilmh^ con- 
straint as 0.130(n/0.98)-i-2± 0.011. While changes of or- 
der 0.08 in tilt are marginally allowed with CMB alone, 
they are strongly disfavored when WMAP is combined 
with the Lyman-alpha forest and galaxy power spectra 
(Seljak et al. 2004). 

Changing the baryon density to Ui,h^ — 0.030 still yields 
a good fit, = 16.2, but increases the inferred ^Imh^ 
to 0.146 ± 0.010. This is not surprising, because higher 
baryon fractions and lower fl„ih^ both increase the ratio 
of large to small-scale power. However, i?o.35 changes only 
to 0.0948 ± 0.0035. This is a 1 cr change in -Ro.35, whereas 
this baryon density change is rejected at 5 cr by WMAP 
(Spergel et al. 2003; Tegmark et al. 2004b). 

As described in §4.2, our model correlation functions 
include two 10% scale-dependent corrections, the first for 
non-linear gravity and the second for scale-dependent bias 
and redshift distortions. Removing the latter changes ^Imh^ 
to 0.148±0.011, a 2 cr change from the baseline. We regard 
ignoring the correction as an extreme alteration. However, 
even this only moves i?o.35 to 0.0985 ± 0.0039. The best 
fit itself is worse, as increases to 19.4. 

The contour plots are based on the covariance matrix 
derived from the mock catalogs. To validate this, we con- 
sider the scatter in the best-fit model parameters among 
the 10 jackknife subsamples (see discussion in § 3.3). The 
jackknifed error in ^Imh^ is 0.011, that in ZJy (0.35) is 4.6%, 
and that in the distance ratio i?o.35 is 3.2%. These val- 
ues are close to those found from the x^-based likelihood 
function (0.011, 4.8%, and 3.7%, respectively). This justi- 
fies the likelihood contours derived from the the covariance 
matrix. It also demonstrates that our results are not being 
driven by one unusual region of the survey. 

4.4. Constraints on dark energy and spatial curvature 

For fixed values of ri„i/i^ and flbh^, the angular scale of 
the CMB acoustic peaks constrains the angular diameter 
distance to z = 1089 to very high accuracy. If one con- 
siders only a simple parameter space of flat cosmologies 
with a cosmological constant, then this distance depends 



only on one parameter, say or f^A (the two must sum 
to unity, and Hq is then fixed by the value of ^Imh^), and 
so the distance measurement constrains firm ^A, and Hq 
to high precision. If one generalizes to larger parameter 
spaces, e.g., adding an unknown dark energy equation of 
state w{z) (Turner & White 1997; Caldwell et al. 1998) or a 
non-zero curvature, then a parameter degeneracy opens in 
the CMB (e.g, Eisenstein, Hu, & Tegmark 1998; Efstathiou 
& Bond 1999). The acoustic scale still provides one high 
quality constraint in this higher dimensional space, but 
the remaining directions are constrained only poorly by 
gravitational lensing (e.g., Seljak 1996; Efstathiou & Bond 
1999) and the Integrated Sachs- Wolfe effect on small and 
large angular scales, respectively. 

With our measurement of the acoustic scale at z = 0.35, 
we can add another high quality constraint, thereby yield- 
ing good measurements on a two-dimensional space (e.g. 
a constant w ^ —1 or a, non-zero curvature). A more 
general w{z) model would of course require additional in- 
put data, e.g. large-scale structure at another redshift, 
supernovae distance measurements, or a Hubble constant 
measurement, etc. 

For the simple space of flat cosmologies with a constant 
w ^ —1, at each value of flmh^ and w, we can find the 
value of ilm (or Hq) that yields the correct angular scale 
of the acoustic peak and then use this value to predict 
Dv{z = 0.35). In Figure 10a, we overlay our constraints 
with the grid of Qm and w infered in this way. One sees 
that flm is well constrained, but w is not. The reason for 
the latter is that flmh^ is not yet known well enough. This 
is illustrated in Figure 11. Were flm.h'^ known to 1%, as is 
expected from the Planck"^*^ mission, then our constraints 
on w would actually be better than 0.1. Figure 10b shows 
the constraints with a grid of Hq and w overlaid. 

Tegmark et al. (2004b) used a Markov chain analysis of 
the WMAP data combined with the SDSS Main sample 
galaxy power spectrum to constrain cosmological parame- 
ters. They found f^^/i^ = 0.145 ±0.014, w = -0.92 ±0.30, 

•^^ http:/ /www.rssd.esa.int/indox.php?project=PLANCK 
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Table 2 



Joint Constraints on Cosmological Parameters including CMB data 



Constant w flat w = —I curved ui = — 1 flat 

Parameter WMAP+Main +LRG WMAP+Main +LRG WMAP+Main +LRG 

w -0.92 ±0.30 -0.80 ±0.18 

■•• ••■ -0.045 ±0.032 -0.010 ±0.009 

Qmh'^ 0.145 ±0.014 0.135 ±0.008 0.134 ±0.012 0.136 ± 0.008 0.146 ±0.009 0.142 ± 0.005 

flm 0.329 ±0.074 0.326 ±0.037 0.431 ± 0.096 0.306 ± 0.027 0.305 ± 0.042 0.298 ± 0.025 

h 0.679 ±0.100 0.648 ±0.045 0.569 ± 0.082 0.669 ± 0.028 0.696 ± 0.033 0.692 ± 0.021 

n 0.984 ±0.033 0.983 ± 0.035 0.964 ± 0.032 0.973 ± 0.030 0.980 ± 0.031 0.963 ± 0.022 



NOTES. — Constraints on cosmological parameters from the Markov chain analysis. The flrst two data columns are for spatially flat models 
with constant w, while the next two are for w = —1 models with spatial curvature. In each case, the other parameters 

and the optical depth t (which we have required to be less than 0.3). A negative means a spherical geometry. The mean values are listed 
with the 1 a errors. The first column in each set gives the constraints from Tegmark et al. (2004b) from combining WMAP and the SDSS 
Main sample. The second column adds our LRG constraints: -Ro.35 = 0.0979 ± 0.036 and and Qmh^ = 0.130(n/0.98)~^'^ ± 0.011. In all cases, 
Qft/i^ is constrained by the CMB to an accuracy well below where we would need to include variations in the LRG analysis. 



= 0.329 ± 0.074, and h = 0.68 ± 0.10 (varying also n, 
flf,, the optical depth r, and a linear bias). Here we use 
the mean and standard deviation rather than the asym- 
metric quantiles in Tegmark et al. (2004b), and we use a 
prior of r < 0.3. Adding the LRG measurement of i?o.35 
and the constraint that fJ^/i^ = 0.130(n/0.98)-i-2±0.011, 
we find = 0.135 ± 0.008, w = -0.80 ± 0.18, f7„ = 

0.326 ± 0.037, and h = 0.648 ± 0.045. We are ignoring the 
small overlap in survey region between SDSS Main and 
the LRG sample. The improvements in w arise primarily 
from the constraint on Vlmh^, while the improvements in 
Q.m come more from the measurement of i?o.35- Table 2 
summarizes these numerical results. 

It is important to remember that constant w models are 
not necessarily good representations of physical models of 
dynamical dark energy and that forcing this parameteri- 
zation can lead to bias (Maor et al. 2004; Bassett et al. 




Fig. 11. — Contours in the space of Q,m, and w. The solid black 
contours show the lines of constant i?o.35 (from 0.090 to 0.106, with 
the central value of 0.098). The dashed red contours show the con- 
tours of constant Qrnh^, using the angular scale of the CMB acoustic 
peaks to set Hq at each (Sim, w) pair. The values of Qmh^ range 
from 0.11 to 0.15, which is the —2 cr to 2 cr range from Figure 7. Un- 
certainties in the value of Cm/i^ significantly impact our constraints 
on w. 



2004). We offer the previous analysis as a means to com- 
pare to the literature, but we prefer our actual distance 
measurements in Table 1 as a model-independent set of 
constraints. 

We next turn to the space of models with two well- 
specified ingredients, namely a cosmological constant (i.e., 
w = —1) and non-zero spatial curvature. The results are 
in Figure 12. Unlike the constraints on w ^ —1, the con- 
straints on the spatial curvature are excellent, of order 1%. 
This is because the distance to z = 1089 is extremely sen- 
sitive to spatial curvature, such that we get excellent per- 
formance by supplying a calibration of the distance scale 
(e.g., with a touch of ^Im) at low redshift. Of course, 
this is in accord with the conventional wisdom that the 
CMB constrains the universe to be nearly flat, but our 
result represents a significant tightening of the angular di- 
ameter distance degeneracy. 

Using the Tegmark et al. (2004b) Markov chain results 
for a w = — 1 cosmology with spatial curvature, the SDSS 
Main sample P{k) and WMAP produces Q.„,h'^ = 0.134 ± 
0.012, VLk = -0.045 ± 0.032, Q.^ = 0.43 ± 0.096, and 
h — 0.57 ± 0.08. Adding the i?o.35 constraint from the 
SDSS LRG results, we find n„,h^ = 0.142 ± 0.011, VIk = 
-0.006±0.011, Vim = 0.309±0.086, and h = 0.679±0.033. 
Adding the further information on Vlmh? drops the val- 
ues to Vtmh'^ = 0.136 ± 0.008, VIk = -0.010 ± 0.009, 
Vim = 0.306 ±0.027, and h = 0.669 ±0.028. Hence, the es- 
sential improvement comes from the measurement of i?o.35; 
with it, we find that the universe is flat to 1% accuracy, 
assuming that w = —1. 

If we require a flat cosmology with w — —1, then the 
Markov chain analysis from the WMAP, Main, and LRG 
data together yields n^h? = 0.142 ± 0.005, 0.^ = 0.296 ± 
0.025, and h = 0.692±0.021. The WMAP data alone is not 
strongly degenerate in this parameter space (Spergel et al. 
2003), although the galaxy data do tighten the constraints 
by roughly a factor of 3. One could also read Vim and 
h directly from the "flat" line in Figure 12. This gives 
Vim = 0.271 ± 0.022 and h = 0.723 ± 0.017 at a fixed 
spectral tilt n = 0.98. The difference occurs because the 
WMAP and SDSS Main pull the value of Vlmh"^ higher 
and because the best-fit tilt in the Markov chain is below 
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Fig. 12. — a) As Figure 10, but now overplotted with models of w = —1 but non-zero curvature. The solid red lines are lines of constant 
curvature, running from —0.03 (closed) to +0.02 (open). The dashed lines are lines of constant f!m, from 0.20 to 0.40. The constraints on the 
curvature of the universe are superb. This happens because we are combining the low-redshift distance scale with the distance to ^ = 1089, 
by which point even tiny amounts of curvature make a big difference, b) As (a), but the dashed lines are now lines of constant Hq. 



n = 0.98. 

4.5. Low-redshift cosmological constraints 

In more general dark energy models, the i?o.35 measure- 
ment will not measure fix or w{z) by itself. However, the 
redshift of the LRG sample is low enough that we can get 
interesting constraints focusing on the path from z = to 
2 = 0.35 rather than z = 0. 35 to z = 1089. We note that 
the combination £)y(0.35)v"m7? has no dependence on 
the Hubble constant Hq, since Dv{0.35) is proportional 
to Hq^ (times a function of all the fi's and w{z)). Fortu- 
itously, this combination is well constrained by our data, 
as these contours lie along the long axis of our constraint 
region. We measure 

/n IP 

A EE £'1/(0.35)-^^-^^ = 0.469 ± 0.017(3.6%). (4) 
0.35c 

This value is robust against changes in the minimum scale 
of data used in the fit (0.471 ± 0.021 for r > 18h-^ Mpc), 
the spectral tilt (0.483 ± 0.018 for n = 0.90), and the 
baryon density (0.468 ± 0.017 for nbh'^ = 0.030). As A is 
independent of a dark energy model, we include its value 
in Table 1. 

If the LRG redshift were closer to 0, then A would sim- 
ply be \/TT^. At z = 0.35, A depends weakly on ^Ik and 
on w{z) over the range < z < 0.35. In detail, for a flat 
universe and constant w, which we denote as wq given the 
low redshift, we have 



A 



^E{zl) 



-1/3 



dz 



E{z) 



2/3 



(5) 



where i;(z) =H{z)/Ho = [n„,{l + zf + ^1^(1 + zf+^'^o] 
and Zi = 0.35. The generalization to curved space-times 
is straightforward. While treating w as a constant for all 
times may be a poor model (Maor et al. 2004; Bassett et 
al. 2004), it is a reasonable approximation for so short an 
interval. In detail, wq is not the value at z = but rather 
some average out to z = 0.35. 

We therefore linearize the expression for A in il„i, fix, 
and Wo to find 

= 0.273 -t- 0.123(1 -t- Wo) + O.UIQk ± 0.025 (6) 



1/2 



This result relies on the acoustic length scale being pre- 
dicted correctly at z ~ 1000 from the CMB measurement 
of rif,/i^ and the LRG measurement of il.mh^, but it is in- 
dependent of the angular acoustic scale in the CMB and 
hence makes no assumption about w{z) at z > 0.35. It will 
depend slightly on unmarginalized parameters such as the 
spectral tilt, the neutrino mass, or other manners of alter- 
ing the LRG value of ^mh^- As we demonstrated in §4.4, 
the CMB acoustic scale is very sensitive to VLk] invoking a 
large \Q,k\ in equation (6) would require large contortions 
in w{z) to maintain the angular location of the acoustic 
peaks. The error in equation (6) is consistent with the 
error on U,m in the constant w Markov chain because the 
uncertainties in wo increase the allowed range of 

5. CONCLUSIONS 

We have presented the large-scale correlation function 
from the SDSS Luminous Red Galaxy sample. This is the 
largest effective volume yet surveyed by a factor of ~ 4 
at small wavenumber. We find clear evidence (3.4 a) for 
the acoustic peak at 100/i~^ Mpc scale. The scale and 
amplitude of this peak is in excellent agreement with the 
prediction from the ACDM interpretations of CMB data 
such as from WMAP. Moreover, the broadband shape of 
the rest of the correlation function gives a measurement of 
the matter density ^Imh^ that matches the CMB findings. 

Before reviewing the quantitative conclusions, we focus 
on the more fundamental ones. The imprint of the acous- 
tic oscillations on the low-redshift clustering of matter is 
a generic prediction of CDM cosmological theory (Peebles 
& Yu 1970; Bond & Efstathiou 1984; Holtzmann 1989; Hu 
& Sugiyama 1996). Our detection confirms two aspects of 
the theory: first, that the oscillations occur at z > 1000, 
and second that they survive the intervening time to be 
detected at low redshift. The small amplitude of the fea- 
tures requires that there exists matter at z 1000 that 
does not interact with the photon-baryon fluid, i.e. dark 
matter. Fully baryonic models or those with extra inter- 
acting matter produce much stronger acoustic signatures 
that would have to be erased by some exotic later process 
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to match low-rcdshift observations. 

In CDM models, large-scale fluctuations have grown since 
z 1000 by gravitational instability. In particular, per- 
turbation theory predicts that small perturbations grow in 
a manner that leaves the Fourier modes of the density field 
uncoupled. This in turn protects the narrow features such 
as the acoustic oscillations as they grow. Non-linear grav- 
itational perturbation theories generically predict mode 
coupling that would wash out the acoustic signature (e.g. 
Fry 1984; Goroff et al. 1986; Jain & Bcrtschinger 1994). 

Hence, the detection of low-redshift acoustic oscillations 
is a stark confirmation of the CDM theory for the growth 
of cosmological structure and the link between the CMB 
anisotropies and the matter perturbations. While the agree- 
ment between recent results on the broadband shape of the 
matter power spectrum (Efstathiou et al. 2002; Tegmark et 
al. 2004a; McDonald et al. 2004), the galaxy three-point 
correlation function (e.g., Feldman et al. 2001), and the 
inferences from the CMB (e.g., Spergel et al. 2003) were 
certainly compelling on this point, we regard the acoustic 
signature as a smoking gun, as its narrowness in real space 
would be diSicult to mimic in alternative models of struc- 
ture formation. This detection confirms the applicability 
of linear cosmological perturbation theory on large scales 
and across a factor of 800 in cosmic expansion. 

The narrowness of the acoustic peak in real space offers 
an opportunity to measure distances to higher redshifts 
(Eisenstein, Hu, & Tegmark 1998; Eiscnstcin 2003; Blake 
& Glazebrook 2003). It is worth noting that this is a cir- 
cumstance where a given improvement in signal-to-noise 
ratio in the clustering statistic makes a super-linear im- 
provement in the distance constraint. One can draw an 
analogy to the determination of the redshift of a galaxy 
with an emission line. A factor of two in signal-to-noise 
ratio can make the difference between detecting the line, 
and hence constraining the redshift to very high precision, 
and not detecting it and having to rely on the spectral 
shape for a low-precision photometric redshift. In the case 
of large-scale clustering, the acoustic scale is not as narrow 
and so the improvement is less dramatic, but we clearly 
benefit at the factor of two level from using the acoustic 
scale rather than the broadband shape of the correlation 
function (i.e. the equality scale). 

In the LRG sample, we measure the acoustic scale to just 
better than 4% precision (1 cr). Comparing this scale to the 
angular scale of the CMB anisotropies gives the distance 
ratio Eo.35 = £'y(0.35)/£)M(1089) = 0.0979 ± 0.0036, 
where £>y(0.35) is defined in equation (2). This distance 
ratio is robust against changes in the broadband clustering 
signal such as via the spectral tilt and against variations in 
our analysis. It is also robust against certain kinds of ex- 
otica, such as adding additional relativistic energy to the 
universe (Eisenstein & White 2004). It does rely on the 
well-understood linear perturbation theory of the recom- 
bination epoch to relate the perturbations in the photons 
to those in the matter. Given this theory, we have mea- 
sured the relative distance between two radically different 
redshifts using a purely geometric method and the same 
physical mechanism. 

This distance ratio is consistent with the familiar cos- 
mological constant cosmology. It is grossly inconsistent 
with the Einstein-de Sitter (Q^ = 1) model, which pre- 



dicts i?o.35 — 0.133 (nominally 10 cr). A model lacking 
dark energy would require fim = 0.70 with fix = 0.30 to 
match the distance ratio. This would require h = 0.90 
and flmh"^ = 0.57 to match the CMB peak location, im- 
plying an age of 8 Gyr. This is in complete disagree- 
ment with the observed shape of the CMB anisotropy 
spectrum, the galaxy correlation function (including these 
LRG data), the cluster baryon fraction (White et al. 1993), 
the observed value of Hq (Freedman et al. 2001), and 
the age of old stars (Krauss & Chaboyer 2003, and ref- 
erences therein), as well as other cosmological measure- 
ments. Hence, our measurement provides geometric evi- 
dence for dark energy. 

The size of the acoustic scale is predicted by very sim- 
ple physics, namely the comoving distance that a sound 
wave can travel between the generation of the perturba- 
tions and the epoch of recombination. In the standard 
cosmological model, this depends only on the matter den- 
sity flmh^ and baryon density fl^h^. The uncertainties on 
flbh"^ from CMB and Big Bang nucleosynthesis are small, 
contributing < 2% to the error on the acoustic scale. How- 
ever, current uncertainties in flmh^ are large enough that 
we need to track the covariance of flmh^ with our distance 
inferences. 

Because the acoustic scale is detected, the LRG data 
alone constrain the equality scale and matter density, i.e. 
we measure flmh^ rather than the more familiar F = flmh. 
With our baseline method and flbh^ = 0.024, wc find 
^mh^ = 0.130(n/0.98)i-2 ± 0.010. This precision is simi- 
lar to that from current CMB measurements (Spergel et 
al. 2003). Importantly, the LRG value agrees with the 
CMB value and with the inference from the clustering of 
the lower redshift SDSS Main galaxy sample (Tegmark 
et al. 2004a), a remarkable cosmological consistency test. 
Because the formal precision of the LRGs is as good as 
the other measurements, we choose to use only the LRGs 
in our fitting. Adding the WMAP information on ilmh^ 
as an external prior would improve the quantitative con- 
straints only slightly, and so we leave it as a cross-check. 
We expect our knowledge of flmh^ to improve rapidly in 
the coming years both from the CMB, with additional 
WMAP data and smaller angle ground-based observations, 
and from large-scale structure, e.g. with improved mod- 
eling of scales below 10h~^ Mpc and the continued data 
collection for the SDSS LRG sample. 

Using the LRG value for il„ih^, we find the distance to 
z = 0.35 to be L»y(0.35) = 1370 ±64 Mpc, a 5% measure- 
ment. Were this at 2; ~ 0, we would have a measurement 
of Hq and Qm, but at z = 0.35 dark energy and curvature 
do matter. The combination _Dy (0.35)-\/Om-ffo/0-35c is 
measured to 4% precision and is independent of Hq . From 
this, we infer = 0.273-|-0.123(1-|-u)o)+0.137^1k±0.025, 
where this wq is the effective value in the range < z < 
0.35. 

Combining with the CMB acoustic scale, we put con- 
straints on more restricted models, either constant w or 
w = — 1 plus curvature. Wc find that our w leverage is 
roughly 0.2. Improvements in knowledge of ilmh^ will 
help significantly. Our leverage on spatial curvature is 
exquisite: we measure Qk = —0.010 ± 0.009. Of course, 
this is a manifestation of the well-known sensitivity of the 
CMB to spatial curvature, but we are breaking the angu- 
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lar diameter distance degeneracy with the best precision 
to date and with a distance ratio that relies on the same 
physics as the CMB. 

It is important to note that because the low-redshift 
acoustic oscillation method measures distances that can 
be quantitatively compared to those from the CMB, the 
method retains sensitivity to phenomena that have more 
effect at higher redshift, such as curvature. Relative dis- 
tance methods, such as supernovae, can only constrain the 
Hubble relation out to the maximum redshift of the sam- 
ple, but absolute methods probe both above and below 
that redshift by using both z = and the CMB as com- 
parison points. 

It is interesting to compare the provenance of the abso- 
lute distance measurements offered by the acoustic scale to 
those of the classical measurements of Hq (e.g. Freedman 
et al. 2001). Once established at 2; > 1000, the acoustic 
scale can be used at low redshift as a standard ruler on 
equal footing to any other. The issue of course is that 
the early Universe (10'^ < z < 10^) is a remote place to 
calibrate one's ruler. There are assumptions about the 
relativistic energy density, the adiabatic nature of the per- 
turbations, the early generation (z > 10"'') of the pertur- 
bations, and the absence of particle decays at 2; < 10^ 
(Eisenstein & White 2004). Many possible alterations cre- 
ate glaring deviations in the CMB anisotropics (e.g. Mood- 
ley et al. 2004). Others, such as small alterations to the 
relativistic density, arc more subtle, at least with present 
data. Our sense is that altering the acoustic scale so as to 
misestimate Hq will require some interesting piece of new 
fundamental physics. The future of direct Hq studies may 
be as a probe of high-redshift particle physics! 

This detection of the acoustic peak at low redshift dra- 
matically confirms several basic assumptions of cosmo- 
logical structure formation theory, but it also points the 
way to a new application of large-scale structure surveys 
for the study of dark energy. Survey volumes of order 
lh~^ Gpc^ offer a reliable standard ruler, whose measure- 
ment across a range of redshifts can determine H{z) and 
Da{z) robustly to percent-level accuracy (Blake & Glaze- 
brook 2003; Hu & Haiman 2003; Sco & Eisenstein 2003). 
Indeed, the available precision improves at higher redshift 
because the acoustic peak is less broadened by non-linear 
structure formation. The observational challenge is to ex- 
ecute these large, wide-field surveys at 2 > 0.5. 
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